Astronomy & Astrophysics manuscript no. maury ' 13492 'text 


©ESO 2010 


January 20, 2010 





o 

G 

o 

(N 



< 

o 

Oh 
6 



On 
cn 

o 
o 



Toward understanding the formation of multiple systems* 

A pilot IRAM-PdBI survey of Class objects 

A. J. Maury 1 , Ph. Andre 1 , P. Hennebelle 2 , F. Motte 1 , D. Stamatellos 3 , M. Bate 4 , A. Belloche 5 , G. Duchene 6 - 7 , and A. 

Whitworth 3 

1 Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/Service d'Astrophysique, C.E. Saclay, Orme des Merisiers, 
91191 Gif-sur-Yvette, France 

2 Laboratoire de radioastronomie, UMR 81 12 du CNRS, Ecole normale superieure et Observatoire de Paris, 24 rue Lhomond, 75231 
Paris, France 

3 School of Physics & Astronomy, Cardiff University, Cardiff, CF24 3 AA, Wales, UK 

4 School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, UK 

5 Max-Planck Institut fur Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

6 Astronomy Department, University of California, Berkeley, CA 94720-3411, USA 

7 Laboratoire d'Astrophysique de Grenoble, Universite Joseph Fourier, BP 53, 38041 Grenoble cedex 9, France 
Received 20 October 2009 / Accepted 13 January 2010 

ABSTRACT 

Context. The formation process of binary stars and multiple systems is poorly understood. The multiplicity rate of Class II pre- 
main-sequence stars and Class I protostars is well documented and known to be high (~ 30% to 50% between ~100 and 4000 AU). 
However, optical / near-infrared observations of Class I/Class II YSOs barely constrain the pristine properties of multiple systems, 
since dynamical evolution can quickly alter these properties during the protostellar phase. 

Aims. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core 
fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at 
the beginning of the main accretion phase. 

Methods. We carried out a pilot high-resolution study of 5 Class objects, including 3 Taurus sources and 2 Perseus sources, using 
the most extended (A) configuration of the IRAM Plateau de Bure Interferometer (PdBI) at 1.3 mm. Our PdBI observations have a 
typical HPBW resolution ~0.3 "— 0.5 " and rms continuum sensitivity — 0. 1 — 1 mjy/beam, which allow us to probe the multiplicity 
of Class protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. 

Results. We detected all 5 primary Class sources in the 1.3 mm dust continuum. A single component associated with the primary 
Class object was detected in the case of the three Taurus sources, while robust evidence of secondary components was found 
toward the two Perseus sources: L1448-C and NGC1333-IR2A. We show that the secondary 1.3 mm continuum component detected 
~ 600 AU south-east of L1448-C, at a position angle close to that of the CO(2-l) jet axis traced by our data, is an outflow feature 
directly associated with the powerful jet driven by L1448-C. The secondary 1.3 mm continuum component detected ~ 1900 AU 
south-east of NGC1333-IR2A may either be a genuine protostellar companion or trace the edge of an outflow cavity. Therefore, our 
PdBI observations revealed only wide (> 1500 AU) protobinary systems and/or outflow-generated features. 

Conclusions. When combined with previous millimeter interferometric observations of Class protostars, our pilot PdBI study 
tentatively suggests that the binary fraction in the ~ 75 - 1000 AU range increases from the Class to the Class I stage. It also 
seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to 
reconcile the low multiplicity rate of Class protostars on small scales with the higher binary fraction observed at later (e.g. Class I) 
evolutionary stages. 
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C3 ' 1- Introduction of cloud core collapse is unclear since it has been argue d that 

most stars may actually form as sing le objects dLadall2006|) . The 

Understanding the formation of multiple systems is a mapr un- ^gamei*. is based on the fact that most low . mass stars (with 

solved problem in star formation research (e.g. | Tohline| | 2002 | ) . ^ < Q 5 M&) are single and that the stellar inMal masfj func . 

While most solar-type (0.5 < M < 2 M ) pre-main sequence tion (IMF) is signi ficantl y more populat ed below 0.5 M than 

(PMS) stars are observed to be in binary systems widitypical aboye ( _ Kroupa 2m . chabrier 20 05). Despite conventional 

separations ~ 10-300 AU (e.g., | Duchene et al. || 2004i |2007l see wisdom ^ it is ther efore conceivable that most low-mass prestel- 

also below), the detailed manner in which individual prestel- lar cores may collapse to single stars 

lar cores fragment (or not) during collapse to form multiple (or It is generally believed that multiple systems form by dynam . 

single) stars is stillnoorly understood and highly debated (see ^ rotationa ii y . dr i V en fragmentation at the end of (or shortly 

i Goodwin etalj | 2007| for a review). Even the typical outcome after) the first coUapse phase of pres tellar cores when the cen- 

tral H2 density reaches « cr „ ~ 3 x 10 10 crrT 3 and the equation 

* Based on observations carried out with the IRAM Plateau de Bure of state of the gas sw itches from isothermality to adiabacity 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG dGoodwin et all 120071) . Purely hydrodynamic SPH simulations 
(Germany), and IGN (Spain). 
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of rotating cloud core collapse show that a very low level of ini- 
tial core turbulence (e.g. E tU rh/E prnv ~ 5%) leads to the forma- 
tion of a multiple system ([Goodwin et alj2004tlHennebelle et all 
l2004t ICommercon et al.ll2008l) . In such SPH simulations, frag- 
mentation is driven by a combination of rotation/turbulence and 
occurs in large ( > 100 AU) disk-l ike structures or "circ umstel- 
lar accretion regions" (CARs - cf. iGoodwin et al.ll2007l) . These 
CARs are not rotationally supported and are highly susceptible 
to spiral instabilities which fragment them into small-Af mul- 
tiple systems with N > 2, typically N ~ 3-4 components at 
radii < 150 AU in the equatorial plane (IGoodwin et al.l 120041: 
iFromang et aT1l2006l) . 

However, a very different outcome is found in simula- 
tions of magnetized core collapse, as shown by recent re- 
sults o btained with MHD c odes using both grid tech- 
niques dFromang et al] 20061: Hennebelle & TevssieJ l2008t 



iMachida et al.l 120051: iMellon & Lil 120081) and the SPH tech- 
nique dPrice & Batd 120071) . These new MHD simulations in- 
dicate that the presence of an even moderate magnetic field 
strongly modifies angular momentum transport during collapse 
and at least partly suppresses core fr agmentation, often l ead- 
ing to the formation of a singl e object. iPrice & Batd d2007l) and 
iHennebelle & Tevssieii d2008l) conclude that binary star forma- 
tion is still possible in the presence of magnetic fields but either 
requires strong initial perturbations or must occur during the sec- 
ond c ollapse phase, after the dissociation of Ho dMachida et al] 
120081) . The systems formed in the latter case a re initially very- 
low-m ass (~ 0.01 M e ), close (~ 1 AU) binaries dBonnell & Batd 
fl99l . which have to grow substantially by accretion during the 
Class 0/Class I phase dBatd l2000|) to match the prop erties of ob- 
served young binary stars (e.g. lDuchene et aO2004l) . Therefore, 
both from an observational and a theoretical point of view, it is 
unclear whether the collapse of a prestellar core typically pro- 
duces one, two, or more stars. 

The multiplicity of solar-type pre-main-sequence stars and 
Class II/Class I young stellar objects (YSOs) is now well doc- 
umen ted and has been investigated at a range of wavelengths 
(e.g.. iGhez et al lll993b ISimon et al.lll995b IPatience et al.ll2002L 
iDuchene et al.ll2004i 120071: IConnellev et al.l l2008h . In particu- 
lar. iPatience et al] d2002TT and iKohler et al] (l2008h showed that 
the binary frequency of (Class II) T Tauri stars ranges from 
~40% to - 60% with a peak in th e separation distribution around 
60+^0 AU Ipuchene et aD d2007l) showed that -32% of Class I 
YSOs have companions in the -50-1000 AU separation range. 
The observed Class I/II binary systems h ave typical mass ra - 
tios q = M second /M main ~ 0.2-1 (e.g. IWoitas et al] l200lh . 
Unfortunately, observations of Class II/Class I YSOs barely con- 
strain the pristine properties of multiple systems since dynami- 
cal evolution c an drastically alter these properties in less than 
~ 10 s yr (cf. iReipurth & Clarkd l200lh . Although still fairly 
uncertain and a matter of debate, the lifetime of the Class 
phase is estimated to be ~ 3 x 10 4 - 10 5 yr dAndre et al] 
l2000t lEvans et alJl2009h. compared to ~ 2 - 5 x 10 s yr for the 
Class I phase dGreene et aljfl994l: lEvans et al]|2009h . Moreover, 
regardless of their precise age and lifetime, Cla ss objects 
are envelope -dominated protostars (M env » M+ - I Andre et al] 
1 1 993l 120001) . while Class I objects tend to have only residual 
protostellar envelopes (M em , < M+ - lAndre & Montmerlell994t 
iMotte & Andrdl2~00lh . Therefore, Class protostars are much 
more likely than Class I objects to retain detailed information 
about the collapse initial conditions and the physics of the bi- 
nary fragmentation process. Furthermore, the circumstellar mass 
reservoir left around Class I sources is generally not sufficient 
to form companions more massive than substellar objects at the 



Class I stage. Probing the multiplicity on scales <1000 AU as 
soon as the Class stage is thus one of the keys to understand- 
ing the bulk of multiple star formation. 

Only subarcsecond mm/submm interferometry can probe 
the inner structure of Class objects, so that little is known 
about their multiplicity on scales < 150 AU. Several inter- 
ferometric studies discover ed a number of wide ( > 1000 AU) 
multiple Class systems ( Loonev et al] 120001 : iBourkd 120011 : 



IChen et al.1 120081: IGirart et al] 1 120091) . but these studies were 



limited by sensitivity to small samples of relatively luminous 
objects. The SMA has recently been used to study several 
Class sources at 1.3 mm and 0.8 mm bu t only at ~ 2" 
(~ 450 AU) resolution dJ0rgensen et al]|2"007h . which does not 
allow the detection of tight multiple systems. Among the eight 
Class sources targeted with the SMA, only one source al- 
ready embedded in a wide separate-envelope system (NGC1333- 
IRAS4A / NGC1333-IRAS4B) was fo und to exhibit a higher 
degree of multiplicity on smaller scales: |j0rgensen et al] d2007h 
showed that IRAS 4A splits into two components separated by 
-2" (450 AU) and IRAS 4B splits into two components sep- 
arated by ~11" (2400 AU). These two new companions were 
detected at both 1.3 mm and 0.8 mm. 

To determine the typical outcome of protostellar collapse and 
constrain binary fragmentation models, high-resolution imaging 
of very young protostars observed as early as possible after the 
end of the first collapse phase are crucially needed. In this pa- 
per, we present the results of a pilot high-resolution survey of 
5 Class objects carried out with the IRAM Plateau de Bure 
Interferometer (PdBI) equipped with new-generation 1 .3 mm re- 
ceivers and using the most extended baselines of the interferom- 
eter (new A configuration). 



2. Observations and data reduction 

2.1. Sample selection 

In order to probe the multiplicity at the Class stage, we con- 
ducted a pilot survey of 5 sources with the Plateau de Bure 
Interferometer (PdBI) in February 2008. The sources were first 
selected based on a distance criterion: they had to be close 
enough so that the PdBI resolution in the most extended con- 
figuration would probe i 150 AU scales, i.e. all sources had 
to be located at d < 250 - 300 pc. The second criterion was 
the source locations in the sky: the sources had to be visible 
from the PdBI in winter and observable with a synthetized beam 
of less than 0.6" (in both directions) with the A array. The se- 
lected sample includes the following five sources: IRAM 04191, 
L1527, L1521F, all located in Taurus at d - 140 pc, and L1448- 
C and NGC 1333-IRAS2A both located in the Perseus cloud at 
d ~ 250 pc (see Table [TJ. These 5 sources have M em , - 0.5- 
4 M , L/,„i < 10 L Q an d are among the y oungest known solar- 



i d are a mong 
dAndre et all 



2000). The main properties 



type Class protostars 
of the 5 Class target sources are summarized below. 



2.1.1. L1448-C 

L1448-C is located in the Perseus molecular clo ud, and was 
first d etected as a 2 cm radio continuum source dCuriel et al.l 
I 990h. associated wi th a strong millimeter continuum source 
( Bachilleretai]|1991l). This is a well-known low-mass Class 
protostar dBarsonv et al]|l998|) drivi ng a powerful, highly colli- 
mated outflow (Bachille r et al.ll 19901) . which has been imaged at 
high resolution (syntheti zed HPBW beam ~2.5"at -90 GHz) in 
CO and SiO with PdBI dGuilloteau et al.lll992t iBachiller et all 
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Source 


(x (JZUUU) 




orK 


Distance 




T 


Ker. 










(PC) 


(Mo) 


(To) 










(1) 


(2) 






(3) 


L1448-C 


03:25:38.87 


+30:44:05.4 


Perseus 


250 


1.6 


5 


(a) 


NGC 1333-IRS2A 


03:28:55.58 


+31:14:37.1 


Perseus 


250 


1.7 


10 


(b) 


IR AM 04191 


04:21:56.91 


+ 15:29:46.1 


Taurus 


140 


0.5-1.5 


0.1 


(c) 


LI 527 


04:39:53.90 


+26:03:10.0 


Taurus 


140 


0.8-1.7 


1.6 


(a), (d) 


L1521-F 


04:28:38.99 


+26:51:35.6 


Taurus 


140 


0.7-4 


0.1 


(e), (f) 



' Star Forming Region with which the Class object is associated 

2) Recent estimat es of the distance t o Perseus range from 220 to 350 pc. Throughout this paper, we adopt a distance of 250 pc for the Perseus 
molecular cloud {Enoch et al. 20Qjj). 

3) References for th e a dopted values of Me „ v a nd L bo[ : (a) iMotte & Andrei ( 120011) . (b) |j0rgensen et alj d2007h . (c) lAndre et all dl999l) . (d) 
lOhashi et al. | jl997l) . (e) lBourke etaD j2006h . (f)|Crapsi et al. (2004]). 



1 1995b . The molecular jet has been extensively studied since its 
discovery, and is known to exhibit very high-velocity features 
(+70 km s 1 - iBachiller et aDfl990lfl995h . 



2.1.2. NGC 1333-IRAS2A 



NGC 1333-IRAS2A was first identified in 450 urn and 850 ^m 



Sandell & Knee 



Rodriguez et al 



continuum observations (ISandell et alJ [1991 
20011) , and is also detected at cm-wavelengths 
1999t iReipurth et alJ l2002h and as a compact 3 mm contin- 
uum source dJ0rgensen et al.l 120041) . CO maps of the IRAS2 
region show two outflows, directed nort h-south and east-west 
(iLiseauet al.ll98 8: Knee & Sandelj2000T) . both originating near 
IRAS2A. Therefore, it has been argued that IRAS2A may be an 
unresolved p rotobinary. However , neither the 2.7 mm BIMA ob- 
servati ons of Loonev et al. d2000h nor the 1 .3 mm SMA observa- 
tions of |j0rgensen et al. d2007l) detect a companion to the source, 
despite a ~3 mJy/beam sensitivity in both cases and beam sizes 
of 0.6" and 2.2", respectively. 



2.1.3. IRAM 04191 

The very low luminosity Class object, IRAM 04191+1522 
(hereafter IRAM 04191), is located in the southern part of the 
Taurus molecular cloud an d was originally d iscovered in the 
millimeter dust continuum dAndre et al]ll999l) . Follow-up ob- 
servations revealed the presence of a CO bipolar outflow and a 
weak 3.6 cm VL A radio cont i nuum source located at its cen- 
ter of symmetry dAndre et al.l 1 19991) . as well as extended in- 
fall and rotation mo tions in a prominent, flattened envelope 
dBelloche et al.l2002h . It is associated with a weak Spitzer source 
and has an estimated accretion luminosity of only L- mt ~ 0.08 L & 
dDunham et alJl2006l) . 



2.1.4. L1527 

L1527 IRS (hereafter L1527), located in the Taurus molecular 
cloud, has been classified as a borderline Class 0/1 object. It is 
obse rved in a nearly ed ge-on configuration (~90° viewing an- 
gle) dOhashi et al.lll997[) and features a large, dense circumstel- 
lar en velope dLadd et aD 199lblChen et al.lll995HMotte & Andrei 
l200l . It also exhibits a prominent bipola r outflow whose lo bes 
are oriented along the east-west direction dParker et al.ll 199 lh . 



2.1.5. L1521-F 

L1521-F, located in the Taurus molecula r cloud, was origi- 
nally classified as an evolved starless core dCodella et allliwk 
lOnishi et al.lll999t ICrapsi et al.ll2004l) . The high central density 
and infall asymmetry seen in the HCO + (3-2) line in dicate an ob- 
ject in the earliest stages of gravitational collapse dOnishi et alJ 
1 1 999h . 12 CO (2-1) observations show no clear evidence of bipo- 
lar outflow emission. Recent Spitzer observations of L1521-F 
detected a low luminosity protostar at mid-infrared wavelengths 
(> 5/an) and in the MIPS 160 ^m data, which led to a reclassifi- 
cation of this objec t as a very low lumin osity Class object (e.g. 
iBourke et al1l2006t iTerebev et ai]|2009t) . Only scattered light is 
detected at near-infrared wavelengths with IRAC, in the form of 
a bipolar nebula oriented east-west which is probably tracing an 
outflow cavity. 



2.2. IRAM Plateau de Bure observations 

Observations of the five sources were carried out at 1.3 mm 
with the IRAM Plateau de Bure Interferometer (PdBI), equipped 
with new-generation receivers in February 2008 (PdBI project 
R068). Broad band continuum emission and 12 CO(2-l) emis- 
sion were observed simultaneously, with the PdBI in its most 
extended configuration (A array with 6 antennas, providing 15 
instantaneous baselines ranging from 24 m to 760 m). The 
proximity of the Taurus and Perseus clouds in the sky allowed 
us to use the same gain calibrators for the two regions, and 
therefore time-share two tracks of ~10 hr on the five sources. 
Each track was divided unequally, depending on the expected 
fluxes of the sources. In particular, a factor of three more time 
was spent integrating on IRAM 04191 and L1521-F than on 
L1448-C. Several nearby phase calibrators (mainly 0415+379 
and 0528+134) were observed to determine the time-dependent 
complex antenna gains. The correlator bandpass was calibrated 
on the strong quasars 3C273 and 3C454.3, while the absolute 
flux density scale was derived from MWC349 and 3C84. The 
absolute flux calibration uncertainty is estimated to be ~15%. 
During the observations, one spectral unit of the correlator was 
tuned to the 12 CO (2-1) line at 230.538 GHz. 

The total bandwidth of this 12 CO spectral unit was 160 MHz, 
with individual channel spacings of 625 kHz (corresponding to a 
velocity resolution of 1.62 km s ). The remaining six windows 
of the correlator were combined to observe the continuum emis- 
sion with a total bandwidth of 1.92 GHz between 229.5 GHz 
and 231.5 GHz. The average system temperature of the 1 mm 
receivers was ~250 K. The typical angular resolution was 0.5"x 
0.3" (HPBW) at the declinations of the targets, while the full 
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Table 2. Rms noise levels and naturally-weighted beam sizes of the final maps 



Source 


230GHz 


230GHz combined (1) 


L,U(z 


1 \ 

-1) 


1U/ OrlZ 




HPBW 


rms 


HPBW 


rms 


HPBW 


rms 


HPBW rms 






(mJy/beam) 




(mJy/beam) 




(mJy/beam) 


(mJy/beam) 


L1448-C 


0.48" x 0.27" 


0.9 


1.68"x 1.39" 


2.8 


0.48" x 0.27" 


10 


4.08"x 3.27" 10 


NGC 1333-IRS2A 


0.49"x 0.32" 


1.16 






0.49"x 0.32" 


20 




IR AM 04191 


0.56"x0.31" 


0.37 


1.37"x0.82" 


0.31 


0.56"x 0.31" 


10 




LI 527 


0.48"x 0.28" 


1.2 


0.88"x 0.79" 


2.0 


0.48" x 0.28" 


20 




L1521-F 


0.49" x 0.27" 


0.12 






0.48" x 0.27" 


8 





(i) L1448-C and L1527: combination of 230 GHz data from project R068 with 218 GHz data from project G080. 



IRAM 04191: combination of 230 GHz data from project R068 with 227 GHz data from Belloche et al. (2002). 



width at half maximum (FWHM) of the PdBI primary beam is 
~ 22" at 230 GHz. 

In this study, we also make use of the 1 .4 mm PdBI observations 
of both L1448-C and L1527, and the 2.8 mm observations of 
L1448-C, carried out in the B, C, D configurations between 
November 1996 and September 1998 (unpublished PdBI project 
G080 by Motte et al.). The typical resolution of these early PdBI 
observations was ~4" at 2.8 mm (107 GHz) and ~2" at 1.4 mm 
(219 GHz). 

In addition, the 2 27 GHz observations obtained by 
iBelloche et al] ( 120021) toward IRAM 04191 in the B, C, D 
configurations of PdBI, which had an HPBW angular resolution 
~2", were also used in combination with our A-configuration 
observations at 230 GHz. 

2.3. Interferomethc data reduction 

All the data were calibrated, mapped, and analyzed with the 
GILDASQ software package. Each map was deconvolved down 
to the th eoretical rm s noise level using the MAPPING CLEAN 
method (IClarkll980h . Natural weighting was applied to the mea- 
sured visibilities, producing synthetized half power beam width 
(HPBW) resolutions ~0.5"x 0.3", as given in Table|2] Note that 
the 3 mm continuum maps of L1448-C and L1527 have sig- 
nificantly larger synthetized beams (HPBW ~ 4"x 3.5") than 
the 1.3 mm maps, as no 3 mm data were obtained with the A- 
configuration. The restored continuum maps have rms values of 
0.12 - 2.8 mJy/beam (see TableO, depending on the integration 
time, array configuration and receivers used during the observa- 
tions. 

Likewise, 12 CO(2-l) data cubes were produced with natural 
uv-weighting, resulting in effective angular resolutions and rms 
noise values reported in the fourth column block of Table [2] 
The 1.3 mm continuum visibilities obtained toward L1448-C 
and LI 527 in the various PdBI configurations were merged to- 
gether in order to produce high spatial dynamic range maps. 
Since the A-configuration data were obtained at a central fre- 
quency of 230.5 GHz, while the BCD-configuration data were 
obtained at a central frequency of 219 GHz, we had to scale the 
BCD-configuration data to 230.5 GHz assuming a spectral index 
a = 3, corresponding to a dust emissivity index/? ~ a-2 ~ 1, for 
the emission of the inner protostellar environment (inner enve- 
lope + disk) traced by our maps. Therefore, the 219 GHz visibil- 
ities were scaled by (230/219) 3 ~1.16 before merging the R068 
and G080 datasets at 230.5 GHz. The respective weights of each 
dataset were adjusted so as to produce the highest possible res- 
olution image, while keeping a low rms noise. The synthetized 

1 Grenoble Image and Line Data Analysis System, 
software provided and actively developed by IRAM 
( http : //www . iram . fr/IRAMFR/GILDAS I 



beam sizes and rms sensitivities of the combined data are given 
in the third column block of Table [2] The resulting maps are 
shown in Fig.|4]and Fig.|7]below, respectively. 
Following the same method, the 227 GHz PdBI visibilities taken 
by Belloche et al. (2002) toward IRAM 04191 and the 230 GHz 
visibilities from our present A-configuration PdBI observations 
were also merged, in order to produce a high spatial dynamic 
range map of IRAM 04191 (see Table© at 1.3 mm. The result- 
ing map is shown in Fig. [7^. 



3. Results of the A-configuration PdBI observations 

3.1. High-resolution 1.3 mm continuum maps 

The 1.3 mm dust continuum maps we obtained with the A-array 
of PdBI toward the 5 sources are shown in Figures 1 and 2. The 
effective spatial resolution of these maps is better than 70 AU 
and 125 AU (HPBW) for the Taurus and Perseus sources, re- 
spectively. All five primary Class targets are detected in these 
high-resolution maps, with signal-to-noise ratios ranging from 
5.4 (IRAM 04191) to 65 (L1448-C) (see Table[3]for flux densi- 
ties and source sizes). 

The three Class sources of Taurus targeted in the present 
study are found to be single in our maps: the main protostellar 
object is the only source detected above the 5<x level in each 
of the maps shown in Fig. 1. A tentative 5<x (~ 6.6 mJy/beam) 
secondary component is detected in the LI 527 map (3.6" west 
of the main source), but it is located at the edge of a dirty lobe, 
which casts doubt on the detection. Moreover, this tentative 
component is not detected in the combined map shown in 
Figure [7] (while such a source should have been detected, given 
the rms noise level of the combined map). Therefore, we will 
not mention this source anymore in the following. 
On the other hand, both Perseus maps (Figs. 2a and 2b) show 
evidence of secondary 1.3 mm continuum components detected 
above the 5cr level. The map of L1448-C reveals a secondary 
1.3 mm continuum source, located ~2.4" south-east from the 
primary L1448-C source (Table [3). It is the only additional 
component detected above the 3<x level in the whole 22" map 
(see Fig. |2] and Fig. [3). The map of NGC 1333-IRAS2A (see 
Fig- EJ a ls° shows a secondary 1.3 mm continuum source (see 
position in Table [3} detected above the 5<x level. The nature of 
these secondary sources is discussed further in Sect. 4 below. 

3.2. n CO(2-1) data 

Compact 12 CO(2-l) emission is detected toward all of the 
targets, except L1521-F. Moreover, significant high-velocity 
12 CO(2-l) emission is detected only toward L1448-C. This is 
due to the fact that the A configuration of PdBI filters out most 
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Table 3. Properties of the 1.3 mm continuum sources detected in the high-resolution PdBI maps 



5 



Source 


Position (1.3 mm) 


Peak flux 


FWHM (I) 


Flux density u) 


Separation {3) 




a (J2000) 


6 (J2000) 


(mJy/beam) 


(arc sec) 


(mjy) 


{ AT T\ 

(AU) 


L1448-C 


03:25:38.87 


30:44:05.3 


59 + 0.9 


0.9 + 0.05 


57 ± 2 




L1448-C / Southl 


03:25:38.95 


30:44:03.2 


11+0.7 


0.7 ±0.1 


11 ± 2 


600 


NGC 1333-IRS2A 


03:28:55.56 


31:14:37.1 


26 ± 1.16 


0.55 ± 0.07 


18 + 2 




NGC 1333-IRS2A / SE 


03:28:56.00 


31:14:33.1 


7 ± 1.16 


0.23 ± 0.04 


4± 1.5 


1900 


IRAM 04191 


04:21:56.90 


15:29:46.1 


2 ± 0.37 


0.28 ± 0.02 


0.9 ± 0.4 




LI 527 


04:39:53.87 


26:03:09.8 


50 ± 1.2 


0.42 ± 0.05 


65 ±3 




L1521-F 


04:28:38.94 


26:51:35.2 


1 + 0.12 


0.65 + 0.009 


0.8 ± 0.2 





<n FWHM diameter is computed from a circular Gaussian fit to the visibilities 

<2) Flux integrated above the 3cr level. 

(3) Projected distance to the primary source. 



of the extended 12 CO(2-l) emission from both protostellar out- 
flows and the parent molecular clouds. 

In the 12 CO(2-l) map of L1448-C, a total of seven compact 
high-velocity CO "bullets" are detected along the bipolar jet 
axis, in both the redshifted and blueshifted lobes (see Fig. |3). 
Two blueshifted bullets are detected above the 3cr level (30 
K.km/s) at LSR velocities ranging from -60 km s _1 to -40 km 
s _1 , and are located 1.44" (360 AU) and 2.69" (670 AU) away 
from the primary source driving the jet, respectively. Five red- 
shifted bullets are detected along the redshifted jet axis, south- 
west of the driving source. These five features have LSR veloc- 
ities ranging from 50 km s _1 to 80 km s , and are located at 
distances from the driving source ranging from 0.9" (125 AU) to 
6.7" (940 AU). The seven high-velocity features detected in the 
L1448-C map are all remarkably well aligned with the axis of the 
molecular jet already mapped at lower resolution by Bachiller 
et al. (1995) with PdBI, and more recently by Jorgensen et 
al. (2007) with SMA. Therefore, we conclude that these high- 
velocity bullets trace the inner part of the jet driven by the pro- 
tostellar source L1448-C. 



4. Nature of the secondary components detected in 
the millimeter continuum maps 

4.1. Sources detected in the vicinity of L1448-C 

4.1 .1 . Secondary 1 .3 mm continuum source 

In the high resolution 1 . 3 mm continuum mapofL1448-C shown 
in Fig. [2^, a secondary source is detected ~2.4" south-east of the 
main source. In the combined 1.3 mm continuum map, this sec- 
ondary component is no longer resolved from the primary source 
(see Fig. |4j>, because its peak flux (11.3 mJy/beam) is only be- 
tween the 3<r and 5<x levels, and it is confused with extended 
1.3 mm emission south-east of L1448-C, which likely arises 
from a cavity in the red-shifted outflow lobe. This secondary 
source is not detected either in the 3 mm map shown in Fig. |4j? 
due to insufficient angular resolution. But it lies close to the 
L1448-C jet axis and is immediately adjacent to the second red- 
shifted high-velocity bullet detected in our 12 CO(2-l) observa- 
tions (only ~0.5" separation - see Fig . [3} . Furthermore, this sec- 
ondary 1.3 mm continuum source coincides wit h the position of 
the fir st SiO(2-l) peak (clump RJ) detected bv iGuilloteau et aO 
(fl99l in the red-shifted outflow lobe. This traces the presence 
of an outflow-induced shock at this position, with a high LSR ve- 
locity offset of +50 km s (see Fig. [5j. Therefore, we conclude 
that the secondary 1.3 mm source detected south-east of L1448- 
C is not a genuine protostellar companion but rather an outflow 



feature directly associated with the high-velocity jet from the 
primary Class object. 



4.1.2. Secondary 3 mm source and adjacent Spitzer source 

Figure|4j5 shows a 3 mm (107 GHz) continuum map of L1448-C 
based on BCD-array data taken with PdBI in 1997 (see Sect. 2.2 
for details). It reveals the presence of a secondary 3 mm con- 
tinuum source, clearly detected above the 5<x level at posi- 
tion (03h25m39.10s, +30°43'57.8"), i.e. -8.4" south-east of the 
main L1448-C source, and which does not have any signifi- 
cant 1 .3 mm counterpart in Figure [4^. A mid-infrared source 
was recently detected with Spitzer 0-6+ 0.2" north of this sec- 
ondary 3 mm source dJ0rgensen et alj|2006l) . and ~7.8" south 
of L1448-C. The small angular separation between the PdBI 
3 mm source and the Spitzer mid-infrared source suggests that 
they are physically related. Neither of our 1.3 mm continuum 
maps shows counterparts to the Spitzer or the 3 mm source above 
the 3cr level, corresponding to upper limits to the 1.3 mm peak 
flux density of 3 mJy/0.48"x 0.27" beam in the high resolu- 
tion 1.3 mm map of L1448-C, and 8.4 mJy/1.68"x 1.39" beam 
in the combined 1.3 mm map. The non detection of the 3 mm 
source at 1.3 mm implies that the spectral index of the emission 
is aHJjJJ™ < 2.0 (after scaling the 1.3 mm and 3 mm fluxes to 
matching beams), excluding the possibility that the whole mil- 
limeter emission is due to dust continuum emission from an em- 
bedded protostellar object (in this case a = J3 + 2, where [3 is 
the dust emissiy ity index, and ~ 3 - 4 is expected - e.g. 

iDent et al.lfl998h . 

If the 3 mm source corresponds to a protostellar object, its 
non-detection in the 1.3 mm continuum emission is unlikely 
to result from interferometric filtering sinc e Class sources 
have strongly centrall y condensed envelopes (lAndre et al.ll2000b 
iMotte & Andril2001h . and are thus expected to be detected in 
a ~1.5" beam. For instance, the 1.3 mm peak flux density of 
L1448-C is ~134 mJy/beam in a ~1.5" beam, so the 3 mm sec- 
ondary source would have to be 17 times weaker at 1.3 mm 
on ~375 AU scales to be undetected (above 3<x level i.e. 8.4 
mJy/beam), which is more than twice the flux ration ~8 mea- 
sured between the two sources computed at 3 mm. To further 
test whether the Spitzer and 3 mm continuum emission may 
originate from a protostellar object, we compared the observa- 
tions with the publicly available grid of model YSO spectral en - 
ergy distributions (SEDs) published bv iRobitaille et al.1 d2007l) . 
We tried to rep r oduce both the mid-infrared fluxes derived by 
|j0rgensen et al.l (120061) for the southern Spitzer source and the 
PdBI 3 mm flux, while keeping the 1.3 mm flux density lower 
than the 5<x detection level (~5 mjy in a 0.37" radius aperture) 
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protostellar object. 
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Fig. 4. (a) Combined 1.3 mm dust continuum map of L1448-C. 
This image was constructed by combining the 230 GHz visibil- 
ities obtained with A-configuration of PdBI with the 219 GHz 
visibilities obtained with the B, C and D-configurations of PdBI. 
The synthetized HPBW is 1.68"x 1.39", and the rms noise is <x 
is ~ 2.8 mJy/beam. The dashed contours are levels of -3<x. The 
full contours are levels of 3<x, 5<x, 8cr and 15<x to 45<r by 15cr. 
The grey circle indicates the position of the secondary 1.3 mm 
continuum source detected in the high resolution map shown in 
Fig. |2 The arrows indicate the direction of the high- velocity jet 
observed in 12 CO(2-l) (see §3.2). (b) 107 GHz continuum map 
of L1448-C. The synthetized HPBW is 4.08"x 3.27", and the 
rms noise is <x ~0.49 mJy/beam. The dashed contour level is - 
2<x. The full contours are levels of 2<x, 5<x, 8<x and 15<x to 85cr by 
15cr. In both panels, the filled ellipses in the lower right corner 
indicate the synthesized HPBW. The star indicates the position 
of the Spitzer source reported bv lj0rgensen et al.l d2006h . 



achieved with PdBI. None of the models explored in this way can 
reproduce the data points properly : the ten best models have to- 
tal^ 2 values ~ 100 - 200 for five data points (compared to total 
X 2 ~ 20-40 with data at the same five wavelengths for L1448- 
C), and show 1.3 mm fluxes which should have been detected 
above the 3<x level in our PdBI observations. The mid-infrared 
emission detected with Spitzer and the adjacent 3 mm emission 
detected with PdBI are therefore unlikely to originate from a 



30"44'00" - 




3 B 25 m 39!0 

RA (J2000) 



38% 



Fig. 5. (a) Image and black contours show the 3 mm continuum 
map of L1448-C (same as Fig. |4). The red and blue contours 
are levels of the SiO(2 -l) line intensity at +65 km s from 
iGuilloteau et al.l (1 19921) . The star indicates the position of the 
Spitzer source reported by Jorgensen et al. (2006), for which no 
counterpart is detected above the 3<x level in the 1.3 mm map 
(see §4.1 for a discussion), (b) Image and black contours show 
the high resolution 1.3 mm map of L1448-C. The red and blue 
contours are levels of SiO(2-l) line intensity at +50 km s . One 
can see that both secondary (1.3 mm and 3 mm) sources coincide 
with peaks of SiO(2-l) emission. 



The 3 mm and Spitzer sources are both located in the walls of 
an outflow cavity (see Jorgensen et al, 2007). They both coin- 
cide with the second SiO(2-l) e mission peak detected toward 
the red-shifted lobe (clump RII in IGuilloteau et al.1l 19921) which 
reveals the presence of an outflow-induced shock at this loca- 
tion, at a high LSR velocity offset of +65 km s _1 (see Fig. [5j. 
Furthermore, these two sources also coincide with a peak in the 
NH3(2,2)/NH3(1,1) ratio, which traces heating due to the inter- 
action between the energe tic outflow and the ambient molecu- 
lar gas dCuriel et al.ll 19991) . We therefore propose that these two 
adjacent sources, loca ted along the je t axis, are in fact shock- 
generated features (cf. iHartig an 2003) associated with the pow- 
erful outflow driven by L1448-C. Compact mid-infrared con- 
tinuum emission along pro tostellar jets has already been ob- 
served: iLefloch et al.1 d2005l) reported the detection of such fea- 
tures along the HH 2 protostellar jet, and argued that the mid- 
IR emission arises from heating of very small grains formed by 
evaporation of dust grain mantles in outflow-induced shocks. 
Furthermore, the lack of a 1.3 mm counterpart to the 3 mm 
source can be explained by the nature of the 3 mm emission. 
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The non-detection of VLA 2 cm emission bv lCuriel et all (1 19901) 
at the position of the 3 mm emission implies a spectral index 
a~ mm > 1.5. Combined with a'- 3mm < 2.0, this suggests that the 

2cm 3mm ' fcte 

3 mm continuum emission is a combi nation of optically thick 
free-free emission (1 < tr|™ < 2 - e.g. lGhavamian & HartigarJ 

1 1998b and optically thin dust continuum emission (a^™ ~ 2 
- 4) associated with a shock in the L1448-C outflow. The fact 
that this outflow-induced shock (traced by the Spitzer and 3 mm 
sources) is not detected at 2 pm can be ascribe d to high visual 
extinc tion toward the southern lobe (Ay ~ 32 - iDionatos et alj 
I2009h . 

Based on this multiwavelength analysis, we conclude that the 
Spitzer mid-infrared source and the PdBI 3 mm source detected 
~8" south of L1448-C do not correspond to protostellar objects 
but are both tracing heating and compression resulting from an 
oblique shock on the outflow cavity wall. 

4.2. Secondary 1.3 mm continuum source detected near 
NGC 1333-IRAS2A 

Our high resolution PdBI observations of NGC 1333-IRAS2A 
allow us to probe the circumstellar environment of this source 
down to ~90 AU scales, with better sensitivity than that achieved 
in previous interferometric observations. While no secondary 
component is detected above the 2<x level within a radius of 
1000 AU (~4") from the primary source, a new secondary 
1.3 mm component is detected 7.7" (~1900 AU) south-east 
of NGC 1333-IRS2A, with a peak flux of ~7 mJy/beam (see 
Fig. |2j5). This secondary 1.3 mm source is located near the 
east-west outflow originating from the vicinity of IRAS2A (see 
Fig. [6}, but is not associated with any known high-velocity CO 
bullet or shock feature. Therefore, the nature of this source 
is unclear: it may either be a genuine protostellar compan- 
ion or an outflow feature. We name this secondary source 
NGC 1333-IRS2A/SE, as it probably lies within the same enve- 
lope as NGC 1333-IRS2A. The weakness of this source explains 
its non-detection in previous millimeter in terferometric studies 
dLoonev et alJl2OQOh[j0r"gensen et alj|2007l) . 

5. Dust continuum emission detected toward the 
primary Class sources 

A detailed discussion of the small scale properties and detailed 
morphology of the dust emission detected toward the primary 
Class sources is beyond the scope of this paper and will be the 
subject of a forthcoming paper. Here, we provide a simple, qual- 
itative description of the slightly extended 1.3 mm continuum 
emission detected with PdBI toward the primary Class sources 
of our sample. 

Our high-resolution 1.3 mm continuum maps exhibit only little 
extended emission. This can be explained by the very high res- 
olution achieved in these interferometer maps and the lack of 
short-spacing data, which filters out most of the extended emis- 
sion from the envelope material. In the combined 1.3 mm and 
3 mm maps, however, the use of multiple array configurations 
allows us to recover some of the extended emission on scales 
ranging from -100 AU to -500 AU. Note that, as the older PdBI 
observations toward L1448-C and L1527 (G080 data taken in 
the B, C, D configurations) are noisier, the combined 1.3 mm 
maps have an improved spatial dynamic range (facilitating, e.g., 
image reconstruction), but higher rms noise values. The com- 
bined 1.3 mm maps of the two sources IRAM 04191 and L1527 
are shown in Fig. [7^ and Fig. [7J3, respectively, while the com- 



SMA map of Jorgensen (2007) 




Fig. 6. SMA map obtained in 12 CO(2-l) by |j0rgensen et al.l 
(120071) toward NGC 1333-IRS2A. The color contours show 
12 CO(2-l) integrated intensity levels in steps of 3<x (see 
lJ0rgensen et al.l d2007l) for details), with the blue contours indi- 
cating emission integrated from -6 to -1 km s relative to the 
systemic velocity, and red contours emission integrated from 
+ 1 to +6 km s _1 relative to the systemic velocity. The two 
lines mark the direction of the protostellar outflows originat- 
ing near IRAS2A. On large scales, the ~ east-west outflow is 
more collimated than the ~north-south bipolar one, and was 
found to be less energ etic than the north-south bipolar outflow 
dKnee & Sandellll2000t) . The two stars mark the position of the 
1.3 mm continuum sources detected in our PdBI observations. 



bined 1.3 mm map of L1448-C was already shown in Fig. [4^ 
above. Details about the combination procedure can be found in 
Sect. 2.3 above. 




Aa (arcsec) Aa (arcsec) 



Fig. 7. Combined 1.3 mm dust continuum maps of the Taurus 
sources IRAM 04191 and L1527. The filled ellipses in the lower 
right corner of the panels indicate the synthesized HPBW beam. 
(a) IRAM 04191. The synthetized beam (HPBW) is 1.37"x 
0.82", and the Icr noise level is ~ 0.31 mJy/beam. The level 
of the dashed contours is -3cr. The solid contours correspond 
to levels of 3cr, 5cr, 7cr, 10cr, and 13cr. (b) L1527. HPBW is 
0.88"x 0.79", and cr is ~ 2.0 mJy/beam. The dashed contour 
level is -3<x. The solid contours are levels of 3cr, 5cr, and 10cr 
to 40cr by 10cr. The four arrows indicate the cavity walls of the 
bipolar outflow driven by L152 7, as traced by the emissi on maps 
of 12 CO(3-2) and HCO + (1-0) dHogerheijde et al.|[l998h . 
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All three combined 1.3 mm continuum maps show extended 
emission features, which correlate well with the outflow cav- 
ity walls delineated by CO line observations of these Class 
objec ts at similar angular resolution (see, e.g., |j0rgensen et al.l 
120071 for a 12 CO(2-l) map of L1448-C with an angular reso- 
lution comparable to that of our combined 1.3 mm continuum 
map). This is particularly clear in the combined 1.3 mm contin- 
uum map of L1527 shown in Fig. |7j5, where the extended dust 
emission around the protostellar source delineates three arms of 
a cross (see also iMotte & Andrei 1200 lh . The cross-like pattern 
seen toward L1527 coincides very well with t he edges of the out- 
flow c avity traced by CO(3-2) o bservations dHogerheijde et al.l 
119981: IChandler & Richerl 120001) . Furthermore, if we compare 
our 1.3 mm continu um map (Fig. Ub) with th e HCO + (1-0) in- 
terferometric map of iHogerheiide et al] d!998l) . we find that the 
cross-like morphology of the dust continuum emission coincides 
very closely with the features detected in HCO + . A similar cross- 
lik e pattern for the du st continuum emission was also observed 
bv lFuller etalldl 995h toward LI 551 -IRS5 . 
This indicates that at least some of the dust continuum emis- 
sion observed in the immediate vicinity of Class protostars is 
caused by column density enhancements due to compression in 
the cavity walls of their outflows. 

6. Discussion: Constraints on the formation of 
multiple systems 

6.1. Multiplicity rate of Class protostars on - 100 AU scales 

All 5 Class protostars observed in the present pilot PdBI survey 
are single on scales between -75 AU and ~1900 AU. The only 
possible companion found is NGC1333-IRAS2A/SE, which is 
located ~1900 AU away from the primary source NGC 1333- 
IRAS2A. We discuss below four possible explanations to the 
non-detection of close protobinary systems in our sample: (1) 
small sample statistics, (2) selection effects in our sample, (3) 
limited mass (and mass ratio) sensitivity, (4) intrinsically small 
multiplicity fraction for Class protostars on ~75-1900 AU 
scales. 

On the first point, even though we observed a small num- 
ber of objects, the probability of drawing five single proto- 
stars from a binary fraction distribution of ~32%, correspond- 
ing to the bin ary fraction of Class I YSOs in the same sep- 
aration range dDuchene et alJ [2007h . is only ~14%. This sug- 
gests that the Class binary fraction may be lower than that 
of Class I and Class II YSOs, but only with very marginal 
statistical significance (1.5cr confidence level) at the present 
stage. I nterestingly, based on a n extensive study of 189 Class I 
sources. IConnellev et al.l (120081) recently speculated that the op- 
posite trend would be observed, namely that the companion star 
fraction should be larger at the Class stage. 
On the second point, the five sources observed in this study were 
selected based on a distance criterion mainly, and belong to two 
different star-forming regions. The sources in our sample have 
bolometric luminosities ranging from 0.1 L to 10 L Q , suggest- 
ing that they span a relatively wide range of final stellar masses. 
Moreover, our sample includes two Class objects embedded in 
clustered environments and belonging to wide systems (L1448- 
C, NGC 1333-IRAS2A), as well as three relatively isolated ob- 
jects (IRAM 04191, L1521-F, and L1527). Albeit limited by its 
small size, our sample thus avoids the most obvious selection bi- 
ases. 

On the third point, our mass sensitivity is directly determined 
by the rms noise achieved in our high-resolution 1.3 mm con- 



tinuum maps. If we assume that the envelope of any puta- 
tive protostellar companion has similar density and tempera- 
ture profiles to the envelope of the primary Class object, 
our observations are sensitive to low circumstellar mass ratios 
(q = M S econd/.Mmam ~ 0.07), except for the faintest two sources, 
IRAM 04191 and L1521-F, toward which we are only sensitive 
to M seC ond/M ma i n ~ 0.6 - 0.9. Assuming Class systems have a 
distribution of circumstellar mass ratios similar to the distribu- 
tion of stellar mass ratio s observed toward T-Tauri binary sys- 
tems dWoitas et alj|200l h. we estimate that the sensitivity of our 
PdBI observations should allow us to detect > 50% of the Class 
binary systems with separations wider than ~100 AU. 
On the fourth and final point, the multiplicity of Class objects 
on scales ~75 - 1000 AU is not well known. While our sample of 
Class protostars does not show any close multiple system (with 
separat ions <1900 AU), the early BIMA 2.7 mm continuum sur- 
vey by lLoonev et al.l d2000l) revealed a higher multiplicity rate 
in their sample at 2.7 mm. Am ong the nine Class objects with 
separate envelopes observed by lLoonev et al.l d2000h . three close 
binary systems were found with separations <1000 AU, lead- 
ing to a binary fraction of ~33% on scales between 100 AU and 
2000 AU. We stress , however, tha t some of the protostellar com- 
panions detected by |Loonev"e t al. (2000) at 2.7 mm could in fact 
be outflow features, like the secondary 3 mm source detected 
in our 3 mm continuum map of L1448-C (Fig. |4j)). One strik- 
ing ex ample is the prototyp ical Class protostar VLA 1623, for 
which lLoonev et alJ d2000t) detected two 2.7 mm components 
separated by 1.11" (i.e ~150 AU), which they interpreted as 
a proto-binary system. Comparing the hig h-resolution 3.6 cm 
VLA image of IBontemps & Andrei d 19971) with the BIMA 2.7 
mm image, it appears that the western 2.7 mm BIMA compo- 
nent may be associated wi th an HH-like object named HH-A by 
IBontemps & Andrei d!997l) . belonging to a series of HH-like cm 
radio continuum knots almost aligned with the outflow axis. This 
suggests that the BIMA 2.7 mm emission observed toward the 
western component (labeled VLA 1623B in lLoonev et alj[2000b 
may be strongly contaminated or even dominated by free-free 
emission, and may not trace the presence of a bona-fide proto- 
stellar companion. In the light of our PdBI findings for L1448-C, 
careful comparison of the results with the location of the out- 
flows from the primary protostars and detection of both 1 mm 
and 3 mm counterparts (so that a spectral index can be derived 
in the millimeter range) are needed before the protostellar nature 
of secondary components detected in millimeter continuum sur- 
veys can be firmly established. 

Combining the re sults of the BIMA 2.7 mm survey of 
lLoonev et al.l d2000l) with our PdBI 1.3 mm results allows us to 
enlarge the sample of Class sources for which the multiplicity 
rate between ~150 AU and ~1000 AU can be discussed. Based 
on the above-mentioned arguments, we do not take into accoun t 
the secondary 3 mm component detected bv lLoonev et al.l d2000h 
close to VLA 1623, as it is probably an outflow feature simi- 
lar to the one observed along the outflow axis of L1448-C. The 
BIMA 2.7 mm observations probe multiplicity on scales > 100 
- 4000 AU for Taurus and Ophiuchus sources, and > 150 - 
5000 AU for Perseus sources. We stress that the enlarged (PdBI 
+ BIMA) sample is not homogeneous because the BIMA and 
PdBI surveys have differing sensitivities, resolutions, and ob- 
serving frequencies. Nevertheless, a simple merging of the two 
samples allows us to draw interesting conclusions. Among the 9 
Class objects mapped with BIMA, two of them are binary sys- 
tems with separations less than 1000 AU: (NGC 1333-IRAS4A1 
/NGC 1333-IRAS4A2) and (IRAS 16293-2422A/IRAS 16293- 
2422B), while no protobinary system is detected at separations 
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less than 500 AU. Since the merged sample (PdBI+BIMA) has 
14 target sources, this leads to estimates of ~14% for the binary 
fraction of Class protostars on ~150 - 1000 AU scales , and 
i 7% for the binary fraction of Class protostars on ~150 - 
550 AU scales. Assuming that the intrinsic binary fraction of 
Class objects in the separation range ~150 - 1000 AU is th e 
same as that of Class I YSOs, i.e. -32% (IConnellev et al.ll2008h . 
the probability of drawing 2 Class binary systems in this sep- 
aration range is only ~9%. Similarly, the probability of draw- 
ing 14 single Class protostars in the separation range 150— 
550 AU, assuming the b inary fraction of Class I YSOs, i.e. ~18% 
(IConnellev et al.ll2008l) . is even lower: ~6%. Therefore, we see 
that combining our sample with that of lLoonev et al.l (l2000h al- 
lows us to strengthen the trend pointed out at the beginning of 
this section, namely that the Class binary fraction may be lower 
than that of Class I YSOs, at least on scales -150 - 550 AU. 
Nevertheless, the trend is only present at the < 1 .9<x confidence 
level in the enlarged sample, and thus remains only marginally 
significant. Clearly, more interferometric observations of Class 
objects taken at comparable angular resolutions would be needed 
to confirm this trend and firmly establish that binary properties 
evolve between the Class and the Class I stage. 

6.2. Comparison with numerical models of binary 
fragmentation 

In this section, we compare the results of our high-resolution 
PdBI observations with the predictions of three published nu- 
merical models of star formation, in terms of multiplicity 
and sp atial structure. The first model dStamatellos & Whitworthl 
l2009h deals with the fragmentation of a massive disk around an 
already formed YS O of compa rable mass. The second hydrody- 
namic simulation dBatd 120091) describes the collapse and frag- 
mentation of a 50 M Q cluster-forming clump, and takes into ac- 
count radiat ive feedback from formed protostellar objects. The 
third model dHennebelle & Fromand2008l:lHennebelle & Ciardl 
12009ft includes the effect of magnetic fields and simulates the 
collapse of an individual cloud core into a protostellar system. 
In order to compare the typical outcomes of these simulations 
with our PdBI observations, the model column density images 
(in g.crrT 2 ) were put to the distance and declination of the Taurus 
and Perseus clouds, then converted into flux density maps (in 
mJy/beam) assuming optical ly thin dust emission at 1.3 mm, 
and /ci.3 mm = 0.01 cm 2 g _1 dOssenkopf & Hennind 1 1994b and 
Td — 10 K for the dust properties. The resulting maps were 
convolved with the typical uv-coverage of the PdBI in A con- 
figuration to produce the synthetic 1.3 mm continuum images 
presented in Figs. 8-10 below. 

6.2.1. Numerical simulations without magnetic fields 

Disk fragmentation model 

Th e hydrodynamic simulations of IStamatellos & Whitworthl 
d2009l) . performed with the SPH code DRAGON, demonstrate 
that the outer parts of massive extended disks are likely to un- 
dergo gravitational fragmentation, thus forming low-mass com- 
panions. Because they fragment rapidly, such massive disks are 
unlikely to be observable beyond the Class phase. We stress 
that protostellar collapse is not modeled in these simulations 
which have no protostellar envelope component. 
Figure [HJa) shows the model column density image resulting 
from an hydrodynamic simulation of a 0.7M o disk around a 
Q.1M G star, as seen ~4500 years after the start of the simulation 



(see Fig. 1 of IStamatellos & Whitwor"thll2009l for more details). 
Figure|8lb) and Figure[8lc) present synthetic 1.3 mm continuum 
images resulting from simulations of A-array PdB observations 
of the model placed at the distance of Taurus and Perseus, re- 
spectively. One can see from the white contour and first black 
contour (corresponding to the average 3cr and 5<x levels achieved 
in our PdBI observations, respectively) that we expect the mas- 
sive, extended disk of the model to be detected as a strong, 
well-resolved structure in IRAM PdBI observations. A circular 
gaussian fit to the visibilities of the synthetic images shown in 
Figures 8b and 8c leads to a FWHM diameter of ~ 3.8" + 0.7"at 
the Taurus distance, and ~ 1.8" ± 0.4"at the Perseus distance: 
this is one order of magnitude greater than the FWHM com- 
puted in the same way for the 5 target sources detected in our 
PdBI maps (see Table 3), which show that all 5 target sources 
are compact when observed with the A array. Furthermore, two 
of the Taurus sources we observed (IRAM 04191 and L1521-F) 
have peak fluxes which are more than one order of magnitude 
weaker than the peak flux in the synthetic image of Fig.[8lb). 

Taken at face value, the refore, our PdBI results are no t 
consistent with the model of Stamatellos & Whi tworthl J2009I 
Note, however, that somewhat less massive (~0.1 M Q ) disks, 
or initially massive disks observed at a later evolutionary stage, 
could be seen as compact structures, at the sensitivity achieved in 
our PdBI observations. One possible explanation for the absence 
of massive extended disks in our observations may be that the 
massi ve disks of the model are short-liv ed (~ 10 4 yr), as pointed 
out by IStamatellos & Whitworthl (120091) . On the other hand, the 
presence of massive, infalling envelopes around the Class ob- 
jec ts we observed (which are not modeled in the simulations 
bv IStamatellos & Whitworthll2009l) should tend to refill any ex- 
tended disk present at early times. A more likely explanation in 
the context of this model is that only ~ 20-30% of all solar- 
type protostars m ay develop massive extended disk s at any time 
in their evolution dStamatellos & Whitworfhir2009l) . Our present 
sample is clearly not large enough to rule out this possibility. 
Therefore, while current numerical simulations of massive, ex- 
tended disks do not satisfactorily reproduce our PdBI observa- 
tions, more numerical simulations and more high-resolution ob- 
servations of Class objects would be needed before robust con- 
clusions can be drawn on this scenario. 



Hydrodynamic model including cloud collapse, disk forma- 
tion, and radiative feedback 

In recent hydrodynamical simulations, iBatd (120091) treats 
both cloud collapse and protostar/disk formation, including the 
effect of radiative feedback from newly formed protostellar ob- 
jects. 

The two model sn apshots presented in Fig. [9^ and Fig. [9}i 
(see also Fig.4 of lBatdl2009l) were processed through the PdB 
simulator to produce the synthetic A-array 1.3 mm continuum 
images shown in Fig. |9j?,e and Fig.[9j;,f, with the models placed 
at the distances of Taurus and Perseus, respectively. One can see 
from the white level and first black level (corresponding to the 
average 3<r and 5cr levels achieved in our PdBI observations, re- 
spectively) that the large disk-like rotating structures produced 
by the models are expected to be detected as strong, extended 
or multiple sources in A-array PdBI observations. A circular 
gaussian fit to the visibilities of the synthetic images shown in 
Figures 9e and 9f leads to a FWHM diameter of ~ 3.4" + 0.7"at 
the Taurus distance, and ~ 1.9" ± 0.5"at the Perseus distance: 
this is again one order of magnitude greater than the FWHM 
diameter measured in the same way for the 5 target sources de- 
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tected in our PdBI maps (see Table 3), which show that all 5 
target sources are compact when observed with the A array. 



We conclude that purely hydrodynamic simulations, even if they 
include radiative feedback which inhibits fragmentation close to 
existing protostellar objects, fail to reproduce our PdBI obser- 
vations, because they tend to form massive extended structures 
(with FWHM ~ 300-500 AU) and/or multiple systems. One 
should bear in mind, however, that the absolute ag e of Class 
objec ts (~ 3 x 10 4 -2 x 10 5 yr) is quite uncertain (cf. lEvans et al.l 
I2009h . which casts some doubt on the time at which the models 
should be compared to our observations. 



6.2.2. Magnetohydrodynamic model 

Adaptive Mesh Refinement (AMR) simulations of cloud 
core collapse and f ra gment ation were carried out by 
iHennebelle & Froman g| d2008l) and Hennebelle & Tevssierl 
d2008l). using t he M HD version of the RAMSES code (see 
iFromang et al.l d2006l) . In these MHD simulations, the for- 
mation of a centrifugally- supported d isk is suppressed 

by magnetic braking (see [Mellon & Lil 120081) . Moreover, 
IHennebelle & Tevssieri d2008l) showed that, for rotating dense 
cores with a magnetic field strength typical of values inferred 
from observations (e.g. Crutcher 1999), fragmentation was 
suppressed by magnetic fields if the initial density perturbations 
were too small. 

Two typical outcomes of these simulations are shown in 
Figure [TUJ Fig. [TUh and Fig.fTOH are model column density im- 
ages obtained in the absence and in the presence of a magnetic 
field, respectively. Fig. [TOb and Fig. [TOfe represent synthetic 
images obtained from these models after convolution with the 
dirty beam of PdBI in A configuration and cleaning, assuming 
a distance d = 140 pc, while Fig. ITDb and Fig. [TOf represent 
the same synthetic images, but assuming a distance d = 250 pc. 
The simulations shown in the upper and lower panels started 
from a rotating, centrally-condensed spherical core 0.013 pc 
in radius, with a density contrast of 10. The total mass of the 
core was 2M Q , and the two simulations started from the same 
initial conditions, except for the initial value of the magnetic 
field (B — in the first case; moderate B-field strength in the 
second case, i.e., initial mass-to-flux ratio fi = 3.2 in units of the 
critical value for collapse). Note that the purely hydrodynamic 
simulation of Fig. [TUh produces a quadruple system within a 
radius of ~ 200 AU, while core fragmentation is completely 
suppressed in the moderately magnetized simulation of Fig.fTOtl. 
which leads to the formation of a single stellar object within 
200 AU. 

The magnetized model shown in the lower panels of Fig. [TOl 
reproduces our PdBI observations of the Taurus Class proto- 
stars quite well: the synthetic image (Fig. ITUfe ) shows a unique 
source, with a peak flux density -3-8 mJy/beam, and a size 
comparable to those observed : a circular gaussian fit to the 
visibilities of the synthetic images shown in Figures lOe leads 
to a FWHM diameter of -0.4" at the Taurus distance (-0.2" at 
the Perseus distance), comparable to the ones computed for 
the sources detected in our PdBI maps, and given in Table 3. 
This comparison suggests that magnetic fields are an essential 
ingredient of the early phases of star formation, as they seem 
to be a plausible way to regulate core/disk fragmentation and 
obtain single objects on scales < 300 AU, similar to what is 
observed in our sample. 



6.2.3. Implications for the formation of multiple systems 

Although our sample clearly needs to be extended before general 
conclusions can be drawn, the non-detection of multiple Class 
systems with separations between -100 A U and 600 AU in 
both our sample and the BIMA sample of lLoonev et al.l (l2000h 
already suggests interesting implications for the formation of 
solar-type multiple systems. 

Two alternative scenarios can be proposed. First, it is possible 
that wide (> 600 AU) multiple systems form during the Class 
phase, and that orbital migration then reduces the separation be- 
tween the protoste l lar components d uring evolution to the Class I 
phase (Bateli 20051 Bate et al.ll2.Q02b . Observations (see, for ex- 
ample IConnellev et al.1 20081) suggest that main-sequence sys- 
tems are tighter than T Tauri systems, and that the latter them- 
selves have separations which are somewhat smalle r than those 
observed toward Class I YSOs dPatience et al.ll2002h . Therefore, 
it is possible that the typical separation of multiple systems de - 
creases in the course of YSO evolution (IConnellev et al.ll2008l) . 
even if such an effect is not yet well understood. The for- 
mation of wi de binary systems at the Class stage has been 
envisa ged by iPrice & Batd (|2007|) and IHennebelle & Tevssieri 
(120081) . Based on their MHD numerical simulations, these au- 
thors showed that, in the presence of large-amplitude initial per- 
turbations, each perturbation develops independently leading to 
the formation of a wide protobinary system, which can then 
gravitationally contract to typical PMS binary separations. 
A second possible scenario would be that multiple systems form 
with tight separations <75 AU, explaining the paucity of sys- 
tems in the - 100-500 AU separation range at the Class stage, 
but that these systems then expand to produce multiple systems 
with wider typical separations at the Class I stage. This alter- 
native scenario is plausible if fragmentation occurs during the 
second collapse phase, a pos sibility which h as b een explored in 
the nu merical simulations of iBonnell d 19941) and lMachida et al.l 
(120071) . One of the conditions for this model to produce binary 
systems with typical separations - 100 - 500 AU is that the ini- 
tially very tight protobinary system must gain sufficient angu- 
lar m omentum by accretion (see e.g. lGoodwin et al1l2004t iBatd 
120001) to increase the separation between the two protostellar 
components by the end of the Class phase. 



7. Summary and conclusions 

We carried out a subarcsecond pilot study of 5 Class objects 
at 1.3 mm with the IRAM Plateau de Bure Interferometer in 
its most extended configuration, in an effort to probe proto- 
stellar multiplicity at separations 50 < a < 5000 AU at the 
beginning of the embedded YSO phase. Continuum emission 
and 12 CO(2-l) emission were observed simultaneously, with 
a typical resolution - 0.3"-0.5"and rms sensitivity - 0.1-1 
mJy/beam, which allowed us to study multiplicity down to 
separations a -50 AU and circumstellar mass ratios q -0.07. 

Our main results and conclusions can be summarized as follows: 

1. All five primary Class protostars (IRAM 04191, L1527, 
L1521-F, L1448-C, and NGC1333-IR2A) were detected in 
the 1.3 mm continuum maps, with signal-to-noise ratios 
ranging from 5.4 to 65.5. 

2. Toward L1448-C, a series of seven high-velocity (vur > 
50 km s~') bullets were detected in 12 CO(2-l), which 
trace the axis of the bipolar jet in both the redshifted and 
blueshifted lobes of the outflow. 



A. J. Maury et al.: Toward understanding the formation of multiple systems 



11 



3. Single 1.3 mm continuum components associated with the ing. The work presented in this paper was stimulated by discus- 
primary Class objects were detected in the case of the three sions held in the context of the Marie Curie Research Training 
Taurus sources, while robust evidence of secondary compo- Network "Constellation" (MRTN-CT2006- 035890). 

nents was found toward the two Perseus sources, L1448-C 
and NGC1333-IR2A. 

4. The L1448-C secondary component lies ~ 600 AU south- 
east of the primary source, at a position angle close to that of 
the CO(2-l) jet axis. We show that it is not a genuine proto- 
stellar companion but rather an outflow feature directly asso- 
ciated with the powerful jet driven by L1448-C. The detec- 
tion of compact millimeter continuum emission originating 
from such an outflow-generated feature emphasizes the need 
to observe outflow-shock diagnostics, before any robust con- 
clusion can be drawn on the nature of secondary components 
detected in the vicinity of protostellar objects. 

5. The nature of the NGC1333-IR2A secondary component, 
detected ~ 1900 AU south-east of the primary source, is as 
yet unclear: it may either be a genuine protostellar compan- 
ion or trace the edge of an outflow cavity. 

6. Altogether, our pilot PdBI survey found only evidence of 
outflow-generated features, and/or wide protobinary sys- 
tems: no multiple system was detected at separations 
a <1900 AU in our sample of 5 Class protostars. 

7. Combinin g our results with th e BIMA survey of 9 Class 
objects bv lLoonev et al.l (l2000l) . we argue that there is no evi- 
dence of multiple protostellar systems in the separation range 
150 < a < 550 AU among an enlarged sample of 14 Class 
protostars. Although the millimeter interferometric observa- 
tions available for this enlarged sample are inhomogeneous, 
they tentatively suggest that the Class binary fraction may 
be lower than that of Class I YSOs, at least for separations 
~150 - 500 AU. This tentative evolution of the binary frac- 
tion from the Class to the Class I stage is present at the 
~1.9cr confidence level in the enlarged sample, and thus re- 
quires confirmation. 

8. Comparison of synthetic model images with our PdBI re- 
sults shows that purely hydrodynamic models of protostellar 
collapse and disk formation have difficulties matching our 
observations, since these models typically produce multi- 
ple components, embedded in large-scale rotating structures, 
which are not observed toward our sample of five Class 
sources. These large-scale rotating structures may be short- 
lived, however, and more observations would be needed to 
draw robust conclusions, given the currently large uncertain- 
ties on the Class lifetime. 

9. Comparison of synthetic model images from magnetohydro- 
dynamic models with our PdBI results shows that magne- 
tized models of protostar formation agree better with our ob- 
servations, as magnetic fields tend to prevent the formation 
of extended disk-like structures and to suppress fragmenta- 
tion into multiple components on small scales (100 AU - 
1000 AU). 

10. However, magnetohydrodynamic models may allow wide (> 
1000 AU) and/or very tight (< 30 AU) multiple systems to 
form during the Class phase. The paucity of multiple Class 
systems with separations 150 AU< a < 600 AU, if con- 
firmed by comparable observations of larger source samples, 
may thus favor binary formation scenarios which involve dy- 
namical evolution of the system separations with time. 
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Fig. 1. High resolution 1.3 mm continuum maps of the Taurus sources. In all panels, the filled ellipse in the bottom right corner 
shows the synthetized HPBW beam, (a) IRAM 04191. The synthetized HPBW is 0.57"x 0.33", and the rms noise is cr ~ 0.37 
mJy/beam. The contour levels are -2cr (dashed), 2cr, Act and 5cr. (b) L1527. The HPBW is 0.48"x 0.28", and cr ~ 1.2 mJy/beam. 
The contour levels are -2cr (dashed), 2cr and 5cr (bold), and 10<x, 30cr, 50cr in white, (c) L1521-F. The HPBW is 0.49"x 0.27", 
and cr ~ 0.1 mJy/beam. The contour levels are -2cr (dashed), 2cr, 5cr (bold), and 8cr, 10cr in white. 
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Fig. 2. High resolution 1.3 mm continuum maps of the Perseus sources. In both panels, the filled ellipse in the bottom right corner 
indicates the synthesized HPBW beam, (a) L1448-C. The HPBW is 0.48"x 0.27", and the rms noise is cr ~ 0.93 mJy/beam. The 
contour levels are -2cr (dashed), 2cr, and 5cr (bold), 10cr to 50cr by 20cr in white, (b) NGC 1333-IRAS2A. The HPBW beam is 
0.57"x 0.33", and cr ~ 1.16 mJy/beam. The contour levels are -2cr(dashed), 2cr, 5cr (bold) ; and 12cr, 20cr in white. The two arrows 
show the positions of the two sources detected above the 5cr level. 
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Fig. 3. The left panel shows the high-resolution 1.3 mm map of L1448-C (same as Fig.[2]but showing a wider area covering most 
of the primary beam). The large dotted circle represents the cleaned area of the map. Note that the large-scale map has not been 
corrected for primary beam attenuation at large distances from the centre (FWHM = 22"). The contour levels are -2cr (light blue 
dashed), 2cr, 8cr and 40cr. The right panel shows a blow-up of the central part of the high-resolution 1.3 mm map of L1448-C in 
greyscale, with superimposed contours of the high-velocity 12 CO(2-l) emission integrated from -60 to -40 km s (blue contours) 
and from +50 to +80 km s (red contours). The first 12 CO(2-l) contour corresponds to the 3cr level (30 mJy/beam); the next 
contour levels are 5<x and 10<x. The double arrow marks the direction of the high- velocity jet detected in 12 CO(2-l). The filled 
ellipse at the bottom of the panels indicates the synthesized HPBW resolution. 
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Fig. 8. (a) Model column density image from the radiative hydrodynamic simulation of Stamatellos & Whitworth (2009) (see text 
for further details), (b) synthetic 1.3 mm continuum image resulting from simulated A-array PdB observations of the model shown 
in (a), assuming a distance d = 140 pc (distance to the Taurus complex), (c) Same as (b), but assuming a distance of d — 250 pc 
(distance to the Perseus complex). The color scale is linear and gives an indicative flux density scale in mJy/beam. The peak flux 
density is found to be ~46 mJy/beam and ~29 mJy/beam in (b) and (c), respectively. In both panels, the white dashed contour 
is showing the typical 3cr detection level (~3 mJy/beam) achieved in our PdB-A observations, while the first plain contour is the 
typical 5cr level. The following black contours are levels of 10, 20 and 40<x. The ellipse in the bottom-left corner represents the 
beam size of PdB-A at the corresponding declination. 
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Fig. 9. Model column density images and s ynthet i c 1.3 mm continuum images resulting from simulated A-array PdB observations 
of two typical outcomes from the model bv lBatd d2009l) (seen 2.6xl0 5 years after the start of the collapse, i.e 1.4 times the initial 
free fall time of the cloud). In each row the first image is the model column density image, while the two images on the right are 
the synthetic images obtained after convolution with the PdBI A-array configuration, (a) is a ~1000 AU wide snapshot, while (d) is 
~600 AU wide, (b) and (e) were produced assuming distances d = 140 pc (similar to the distance of the Taurus complex), (c) and 
(f) were produced assuming d = 250 pc (similar to the distance of the L1448 complex). In all panels, the contours are levels of 3, 
5, 10, 20 and 40<x, as achieved in our observations. The color scale is linear and gives an indicative flux density scale in mJy/beam. 
The dashed contour is showing the typical 3<x detection level (~3 mJy/beam) achieved in our PdB-A observations. The ellipse in 
the bottom-left corner represents the beam size of PdB-A at the corresponding declination. 
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Fig. 10. Synthetic 1.3 mm continuum images result ing from simulated A-array Pd B observations of two typical outcomes from the 
simulations of magnetized core collapse from e.g. iHennebelle & Tevssierl (120081) (see text for further details). Upper row: Panel 
(a) shows a model snapshot view of the column density distribution of the inner part of a protostellar system obtained ~ 10 4 yr 
after the beginning of collapse in purely hydrodynamic (B = 0) simulations of cloud core collapse. Panel (b) shows the synthetic 
PdB-A 1.3 mm dust continuum image produced from the purely hydrodynamic model (a), put to the distance of the Taurus complex 
(d = 140 pc). Panel (c) shows a synthetic image similar to (b), but put to the distance of the Perseus complex (d = 250 pc). Lower 
row: Panel (d) shows a snapshot view of the column density distribution of the inner part of a protostellar system obtained from 
MHD simulations starting from the same initial conditions as in panel (a), except for a non-zero magnetic field (whose initial value 
is 1/2 of the critical field strength required to prevent collapse). Panel (e) shows the synthetic PdB-A 1.3 mm dust continuum 
image produced from the MHD model shown in (d), and put to the distance of the Taurus complex (d = 140 pc). Panel (f) shows 
a synthetic image similar to (e), but put to the distance of the Perseus complex (d = 250 pc). In all of the synthetic images, the 
color scale is linear and gives an indicative flux density scale in mJy/beam, assuming Tj — 10 K and /ci.3„„„ = 0.01 cm 2 g _1 for the 
dust properties. The white dashed contour shows the typical 3<x detection level (~3 mJy/beam) achieved in our PdB-A observations, 
while the following black contours are: 5<x (thick), 10cr, 20cr and 40cr. 
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